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Abstract 
The flux of muons induced by atmospheric neutrinos has been measured with the 
qj ■ MACRO detector. Different event topologies have been detected, due to neutrino 

interactions in the apparatus and in the rock below it. The upward-throughgoing 
muon sample is the larger one and is generated by neutrinos with a peak energy 
of ~ 100 GeV. The observed upward-throughgoing muons are 26 % fewer than 
expected and the zenith angle distribution does not fit with the expected one. As- 
suming v^ <-> v T neutrino oscillation, the angular shape and the normalization 
factor suggest maximal mixing and Am 2 of a few times 10~ 3 eV 2 . Also the other 
event categories induced by low-energy neutrinos (peak energy ~ 4 GeV) show a 
deficit of observed events. The value of this deficit and its uniformity with respect 
to the angular bins are in agreement with the oscillation parameters suggested by 
the analysis of the upward-throughgoing muon sample. 



Les Rencontres de Physique de la Vallee d'Aoste 
February 28 - March 6, 1999, La Thuile, Italy 



1 Introduction 

The interest in precise measurements of the flux of neutrinos produced in cosmic ray 
cascades in the atmosphere has been growing over the last years due to the anomaly 
in the ratio of contained muon neutrino to electron neutrino events K The anomaly 
finds explanation in the scenario of neutrino oscillation. The effects of neutrino 
oscillations appear also at higher energies, as reported by many experiments A 
The flux of muon neutrinos in the energy region from a few GeV up to hundreds of 
GeV is inferred from measurements of upgoing muons. The flux of upgoing muons 
is reduced as a consequence of v oscillations. A clearer signature of v oscillations 
in the range 10~ 3 < Am 2 < 10~ 2 eV 2 is connected with the dependence of the 
reduction on the polar angle 9 with respect to the zenith. The reduction in the 
number of events is stronger near the vertical than near the horizontal directions 
due to the longer pathlength of neutrinos from production to observation near the 
nadir (Fig. [J]). 
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Figure 1: Path of atmospheric neutrinos versus cos 9. The upward going v's (neg- 
ative values of cos 9) travel for thousands of kilometres before reaching an under- 
ground detector. 

Furthermore the flux of atmospheric muon neutrinos in the region of a 
few GeV can be studied looking at muons produced inside the detector and muons 
externally produced and stopping inside it. If the atmospheric neutrino anomalies 
are the results of neutrino oscillations, it is expected a reduction in the flux of 
upward-going low-energy atmospheric neutrinos of about a factor of two, but without 
any distortion in the shape of the angular distribution. 

In this paper the MACRO flux measurements are updated for high and 
low energy muon neutrinos. The new results are in agreement with the previous 



ones 



lil' , and the indication for the v oscillation hypothesis is now stronger. 



2 The MACRO detector 



The MACRO detector \=-^) is located in the Gran Sasso Laboratory, with a minimum 
rock overburden of 2700 hg/cm 2 . It is a large rectangular box (76.6 x 12 x 9.3 m 3 ) 
divided longitudinally in 6 supermodules and vertically in a lower and an upper 
part, called attico. The active elements (see Fig. |]) are liquid scintillation counters 
for time measurement and streamer chambers for tracking, with 27° stereo strip 
readouts. In the lower half of the detector streamer tube planes are alternated 
with trays of crushed rock absorbers. The attico is hollow and it is used as a work 
area in which the electronics racks are placed. The streamer tube system allows to 
achieve a tracking resolution for muons in the range 0.2° - 1° as a function of the 
track length. Hence the angular uncertainty is lower than the angular spread due 
to multiple scattering in the rock for a muon coming from above. The scintillator 
system consists of horizontal and vertical layers. The scintillator resolutions for 
muons are about 0.5 ns in time and 11 cm in position. 

Thanks to its large area, fine tracking granularity and electronics symmetry 
with respect to upgoing and downgoing flight directions, the MACRO detector is 
a proper tool for the study of upward-travelling muons, generated by external v 
interactions. Its mass (~ 5.3 kton) allows also to collect a statistically significant 
sample of neutrino events due to internal interactions. 

3 Neutrino events in MACRO 

Fig. |2] displays the different kinds of neutrino events here analyzed. Most of the 
detected particles are muons generated in v^ Charged Current interactions. In 
Fig. H the parent neutrino energy distributions for the different event topologies 
are shown : 

• Up Through - These tracks come from interactions in the rock below MACRO 
and cross the whole detector (E^ > 1 GeV). The time information provided by 
scintillator counters allows to know the flight direction by means of the time- 
of-flight (T.o.F. ) method. The data have been collected in three periods, 
with different detector configurations. In the first two periods (March 1989 - 
November 1991, December 1992 - June 1993) only lower parts of MACRO 
were working B\ In the last period (April 1994 - February 1999) also the 
attico was in acquisition. 



• In Up - These partially contained events come from v interactions inside the 
apparatus. Also in this case the analysis algorithm is based on the T.o.F. 
calculation made possible by the attico scintillator layers. Hence only the data 
collected with the attico (live-time ~ 4.1 years) are useful for this analysis. A 
not negligible percentage (13 %) of events is expected to be induced by Neutral 
Currents or v e CC interactions. 

• Up Stop + In Down - This sample is composed by two subsamples : ex- 
ternal interactions with upward-going track stopping in MACRO (Up Stop), 
interactions in the lower part of the detector with a downgoing track (In 
Down). These events are recognized by means of topological criteria and the 
lack of time information makes difficult to distinguish the two subsamples. Up 
to now we consider them as an unique sample. Assuming that neutrinos do 
not oscillate, it is expected that the number of Up Stop is almost equal to 
the number of In Down and the contribution of Neutral Currents and v e CC 
interactions is ~ 10 %. The attico is used also in the analysis of this sample as 
a veto for downgoing tracks. Therefore the analyzed data are those collected 
with the whole detector with the same effective live-time of the In Up sample. 

4 Analysis procedure 

The T.o.F. method uses the formula 

1 cxj^-Ti) 

PL' { ' 

where T% and T 2 are the times measured in lower and higher scintillator counters, 
respectively, and L is the path between the two counters. Therefore 1/(3 results 
roughly +1 for downgoing tracks and —1 for upgoing tracks. In the sample of 
throughgoing muons it happens that a track hits 3 scintillator layers. In this case 
the time measurements are redundant and 1/(3 is calculated by means of a linear fit 
of the time as a function of the pathlength. 

Several cuts are imposed to remove backgrounds from radioactivity and 
showering events which may cause failure in time reconstruction. Another back- 
ground is connected to photonuclear processes such as fj,N — > fiirX where low-energy 
upgoing particles are produced at large angles by downgoing muons \±H) . The re- 
quirement of a minimum range of 200 g/cm 2 in the apparatus is applied to the 
Up Through sample in order to reduce drastically these low-energy tracks which 
mimic neutrino induced events when the downgoing muon misses MACRO. 
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Figure 2: Cross section view of the detector and topology of the neutrino induced 
events. The stars, the dashed lines and the full lines indicate scintillator hits, neu- 
trino paths and charged particle paths, respectively. 



After all analysis cuts the signal peaks with 1//3 ~ —1 are well isolated for 
the first two samples (see Fig. |j). 

The Up Stop + In Down events are identified via topological constraints. 
The main requirement is the presence of a reconstructed track crossing the bottom 
scintillator layer. All the track hits must be at least 1 m far from the supermodule 
walls. The criteria used to verify that the event vertex (or /i stop point) is inside 
the detector are similar to those used for the In Up search. The probability that an 
atmospheric muon produces a background event is negligible. To reject ambiguous 
and/or wrongly tracked events which survive automated analysis cuts, real events are 
directly scanned with the MACRO Event Display. Because of the low energy of these 
events the minimum range of 200 g/cm 2 in the apparatus is not required. Therefore 
the background due to upward going charged pions is estimated via simulation and 
subtracted on a statistical basis. 

Expected rates and angular distributions have been estimated assuming 

UD HI 



the atmospheric v flux calculated by the Bartol University group 
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Figure 3: Distributions of the parent neutrino energy giving rise to different kinds of 
events, estimated by means of Monte Carlo simulation using the same cuts applied to 
the data. The average energy is ~ 100 GeV for Up Through sample and ~ 4 GeV 
for In Up and Up Stop + In Down samples. 

timate of v cross-section was based on GRV94 LkaJ parton distribution set, which 
varies by +1 % the prediction with respect to the Morfin and Tung li^ parton 
distribution used in the past by the MACRO Collaboration. For In Up and 



Up Stop + In Down samples also low-energy effects Lk3/ have been taken into 
account. The propagation of muons in the rock was taken from Lohmann et aZLLii/ . 
The uncertainty on the expected muon flux is estimated 17 % for Up Through 
events and 25 % for the other events. The apparatus and the data acquisition are 
fully reproduced in a GEANT \=2v based Monte Carlo program. Real and simulated 
data are analyzed by means of the same procedure. Particular care has been taken to 
minimize the systematic uncertainty in the detector acceptance simulation. For the 
Up Through sample different acceptance calculations, including separate electronic 
and data acquisition systems, have been compared. For each sample, two different 
analyses have been performed getting the same results. Furthermore trigger and 
streamer tube efficiency, background subtraction, effects of analysis cuts have been 
in detail studied. The efficiency of the visual scanning for the Up Stop + In Down 
sample has been estimated by analyzing real and simulated events after a random 
merging. The systematic error on the total number of events due to the acceptance 
has been estimated 6 % for Up Through sample. The uncertainty is higher (10 %) 
for low-energy samples because it depends strongly on data taking conditions, anal- 
ysis algorithm efficiency and detector mass. 
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Figure 4: 1//3 distributions after all analysis cuts for throughgoing (first plot) and 
partially contained (second plot) tracks. The data are those collected when also the 
attico was in operation. The neutrino-induced signals are close to 1/(3 = — 1. The 
shaded part in the first plot indicates events whose T.o.F. is calculated by means of 
a fit of three times. In the second plot the events with 1/(3 > are downgoing 
atmospheric muons stopping in the detector. 

5 High energy sample - Analysis results 



In the Up Through sample 642 events are in the signal range (—1.25 < 1/(3 < 
—0.75). Looking at the distribution of events outside the signal peak, we estimate 
a background contamination of 12.5 ± 6 events. Furthermore 10.5 ± 4 events are 
expected to be upgoing charged particles produced by downgoing muons. Finally we 
expect that 12 ±4 are not true Up Through events because they are generated by 
neutrino interactions in the very bottom layer of MACRO scintillators. Hence, re- 
moving these backgrounds, the observed number of Up Through muons integrated 
over all zenith angles results 607. 

For this sample 825 events are expected and the ratio of observed to ex- 
pected events is reported in Table [J]. In Fig. [5] the cos 6 distribution of the measured 
flux is shown compared with the expected ones assuming stable or oscillating neu- 
trinos. The data error bars are the statistical errors with an extension due to the 
systematic errors, added in quadrature. The observed zenith distribution does not 
fit well with the hypothesis of no oscillation, giving a maximum y 2 probability^ of 
only 0.35 %. Combining normalization and angular shape the probability is still 



x The data are normalized to the predicition. The last bin near the horizontal is not taken into 
account because of the higher acceptance uncertainty and the background connected to scattering 
of quasi horizontal downgoing muons. 
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Figure 5: Measured and expected fluxes are shown for the Up Through sample with 
a muon energy threshold of 1 GeV. The solid curve and the shaded region show the 
expectation for no oscillations and its uncertainty (the uncertainty on the predicted 
shape is much smaller). The dashed line shows the predicition assuming v^ <-> u T 
oscillation with maximal mixing and Am 2 = 0.0025 eV 2 . 



very low (0.36 %). 

These results can be explained in the scenario of v oscillation. Assuming 
v^ <-> v T the lower x 2 value for the angular distribution in the physical regionF] is 
12.5 with maximal mixing and Am 2 ~ 0.0025 eV 2 . In the first plot of Fig. ^ the 
independent probabilities for normalization and angular shape and the combined 
probability are shown as functions of Am 2 assuming maximal mixing. It is notable 
that total number of events and cos 9 distribution indicate very close values of 
Am 2 . The second plot of Fig. |^ shows the same probabilities for v^ ^-> v sterile 
oscillation \=^ tM) . The maximum of the combined probability is 36.6 % for v^ <-> v T 
and 8.4 % for oscillation into sterile neutrino. 

Fig. |7| shows the confidence regions at the 90 % and 99 % confidence levels 
based on application of the Monte Carlo prescription by Feldman et al 

2 The best x 2 value is 10.6 in the unphysical range (sin 2 28 mix = 1.5). 
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Figure 6: Probabilities of Up Through muon results assuming v oscillation with 
maximal mixing. The dotted line is referred to the total number of events, the dashed 
line to the \ 2 °f the normalized angular distribution and the continuos line to their 
combination. The v^ ^> v T oscillation is assumed for the first plot, the v^ <-» v s teriie 
for the second one. 

sensitivity of the experiment is plotted. The sensitivity is the 90 % contour which 
would result when the data are equal to the Monte Carlo prediction at the best-fit 
point. 



6 Low energy samples - Analysis results 

In the In Up sample the uncorrelated background is estimated from the 1//3 distri- 
bution. After the background subtraction 116 events are accepted. The prediction 
is 202 In Up events. In the Up Stop + In Down sample 193 events survive the 
analyis cuts and the visual scanning while 274 events are expected. 

The ratios of the observed number of events to the prediction and the 
angular distributions of both samples are reported in Table |I] and in Fig. |8|. The 
low-energy v^ samples show an uniform deficit of the measured number of events 
over the whole angular distribution with respect to the predictions based on the 
absence of neutrino oscillation. 

We note the agreement between the results for low-energy and Up Through 
events. Assuming the oscillation parameters suggested by higher energy sample, it is 
expected a ~ 50 % disappearance of v^ in In Up and Up Stop samples because of 
the neutrino path (thousands of kilometres) . No flux reduction is instead expected 
for In Down events whose neutrino path is of the order of tens of kilometres. The 




Figure 7: The results of the combined (normalization + angular shape) analysis for 
Up Through sample are shown in terms of confidence regions at the 90 % and 99 % 
levels assuming v^ <-> v T . Since the best probability is outside the physical region 
the confidence interval regions are smaller than those expected from the sensitivity 
of the experiment. 
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ratios and the angular distributions estimated assuming the v oscillation are also 
reported in Table [l] and in Fig. || 

In order to reduce the effects of uncertainties coming from neutrino flux 
and cross section, the double ratio 



(2) 



has been studied. A residual theoretical error (5 %) survives, due to small differences 
between the energy spectra of the two samples. Because of some cancellations the 
systematic uncertainty is also reduced to ~ 6 %. The value of the double ratio over 
the zenith angle distribution is shown in Fig. [| Assuming the oscillation parameters 
(sin 2 26^3; = 1, Am 2 = 0.0025 eV 2 ) suggested by Up Through sample, the points 
are compatible with 1. 

The probability to obtain a sum double ratio at least so far from 1 is ~ 6% 
assuming no oscillation and the Bartol flux as the parent v flux. 



Topology 


Ratio 


Statist. 
error 


Syst. 
error 


Theor. 
error 


Ratio with 
v oscillation 


Up Through 

In Up 

Up Stop + In Down 


0.736 

0.57 
0.71 


0.031 
0.05 
0.05 


0.044 
0.06 
0.07 


0.12 
0.14 
0.18 


1.058 

0.98 

0.92 



Table 1: Ratios of observed on expected number of events for different event topolo- 
gies. In the last column the ratio is calculated assuming v^ <-* v T oscillation 
with the parameters suggested by the Up Through angular shape (sin 2 29 m i X = 1, 
Am 2 = 0.0025 eV 2 ). 



7 Conclusions 

The flux and the shape of the zenith distribution for the Up Through sample 
favour i/p <-► v T oscillations. The experimental data have a 36.6 % probability 
assuming oscillation against 0.36 % assuming stable neutrino. 

Therefore the new data confirm the MACRO published results B* \l$ with 
increased probability for the oscillation hypothesis. 

The low-energy neutrino events are fewer than expected and the deficit is 
quite uniform over the whole angular range. Also these results suggest oscillation 
with maximal mixing and Am 2 of a few times 10 -3 eV 2 . 

The combined analysis of high and low energy data is in progress. Presently 
we stress the strong coherence of the MACRO results in different energy ranges and 
with different event topologies. 
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